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Abstract 

We propose an effective method to optimize the working parameters (WPs) of microwave-driven quantum logical gates imple- 
mented with multi-level physical qubits. We show that by treating transitions between each pair of levels independently, intrinsic 
gate errors due primarily to population leakage to undesired states can be estimated accurately from spectroscopic properties 
of the qubits and minimized by choosing appropriate WPs. The validity and efficiency of the approach are demonstrated by 
applying it to optimize the WPs of two coupled rf SQUID flux qubits for controlled-NOT (CNOT) operation. The result of 
this independent transition approximation (ITA) is in good agreement with that of dynamic method (DM). Furthermore, the 
ratio of the speed of ITA to that of DM scales exponentially as 2 n when the number of qubits n increases. 

PACS numbers: 03.67.Lx, 85.25. Dq, 89.70.+C 
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A practical quantum computer would be comprised of 
a large number of coupled qubits and the coupled qubits 
must be kept in high-degree quantum coherence states 
for sufficiently long time. During the past decade, sig- 
nificant progress has been made on physical implemen- 
tation of quantum computation. High-degree quantum 
coherence has been demonstrated experimentally in sys- 
tems such as trapped ions H, 0, nuclear spins 0, Q, 
atoms in optical resonators [a], and photons in microwave 
cavities p]. However, it seems quite difficult to real- 
ize a large number of coupled qubits using these sys- 
tems. Meanwhile, solid-state qubits are of particular 
interest because of their advantages of large-scale inte- 
gration, flexibility in design, and easy connection to con- 
ventional electronic circuits [jj. Of those, qubits based 
on superconducting devices have recently attracted much 
attention || as manipulation of quantum coherent states 
being successfull y d emonstrated in a variety of single 
qubits Ifl EH Ell Il2t EH | and coupled two-qubit systems 
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However, the solid-state qubits demonstrated in exper- 
iments so far all have relatively short coherence time and 
high probability of gate errors H El El El El El El 
Il7j . One of the causes of these problems is extrinsic gate 
error arising from interaction between the environment 
and qubits resulting in decoherence, such as dephasing 
and relaxation |8}. Another cause is intrinsic gate error 
resulting from population leakage to undesired states due 
to the typical multi-level structures of solid-state qubits 
[lii. Eel ] . The intrinsic gate error is crucial since it not 
only contributes to additional decoherence but also de- 
termines the ultimate performance of the quantum gates 
and cannot be eliminated by reducing the environment 
caused decoherence. 

The leakage of a gate can be characterized quantita- 
tively by summing up the maximum transition probabil- 
ities to all undesired states of a multi- level qubit |lfltl2f|. 
It depends strongly on energy level structure and transi- 
tion matrix elements, i.e., spectroscopic properties of the 



multi-level qubit determined completely by device pa- 
rameters (DPs) and external control parameters of the 
qubit which we call working parameters (WPs) for sim- 
plicity. For instance, inductance (capacitance) of the su- 
perconducting flux (charge) qubit is a DP while external 
flux (gate voltage) is a WP. For the multi-level qubit with 
given DPs, the gate leakage is sensitive to their WPs and 
thus can be minimized by the use of appropriate WPs. 

Conventionally, the leakage is calculated from tran- 
sition probabilities to all undesired states by numeri- 
cally sol ying th e time-dependent Schrodinger equation 
(TDSE) 20, 21] . However, this kind of dynamic method 
(DM) not only needs complicated numerical algorithms 
but also a large amount of computational resources. For 
instance, optimizing the WPs of a many-qubit network 
needed for a practical quantum computer using DM may 
only be possible with ad hoc powerful quantum comput- 
ers in the future. Thus a much faster approach is highly 
desirable. 

In this Letter, we propose a very fast method to min- 
imize the leakage of microwave-driven quantum logical 
gates by choosing appropriate WPs in a system of multi- 
level qubits with their DPs given in prior. We consider 
the case of weak microwave fields only since strong fields 
usually cause many additional types of intrins ic g ate er- 
rors and thus should be avoided in general |2(j. Our 
method is based on an independent transition approxi- 
mation (ITA) in which transitions in multi-level qubits 
are treated independently. The leakage is estimated us- 
ing the spectroscopic properties obtained by solving the 
eigenvalue equation of and minimized by optimizing the 
WPs of the qubits. The method is applied to minimize 
leakage of controlled-NOT (CNOT) gate implemented 
with coupled rf superconducting quantum interference 
device (SQUID) flux qubits. The result is in good agree- 
ment with that obtained from DM. More importantly, 
the ITA is scalable as the number of qubit n increases 
because the ratio of the speed of ITA to that of DM 
scales exponentially as 2™. 
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A microwave-driven gate is realized via coherent tran- 
sitions between the computational states of qubits inter- 
acting with microwave fields. In general, correlation and 
interference between transitions cannot be ignored and 
transition probability can only be computed accurately 
by solving TDSE. However, in weak fields, which is the 
case considered here, only the transitions between levels 
with which the microwave field is resonant or nearly reso- 
nant are significant. For a given level, if the level spacings 
between it and all other levels are sufficiently different, 
the correlation and interference between the transitions 
from this level have negligible effect. Hence, each tran- 
sition is expected to take place independently and the 
two levels involved can thus be treated as if they are 
isolated from the others. For each of the two-level sub- 
systems interacting with a rectangular pulse, the transi- 
tion probability can be approximated by an analytical ex- 
pression using the rotating- wave approximation (RWA). 
Assuming the rectangular pulse is e(t) = eo cos (cot), 
the interaction between the qubit and the pulse is then 
V(t) = —fj,-e(t), where eo and 10 are the amplitude and 
frequency of the pulse and fi is the dipole moment oper- 
ator of the qubit. The maximum transition probability 
from state i to j in an A-photon process, Vij, is 22] 

P« = fi?-/(P£ + «?•)> (1) 

where, = 2/iy-eoA Jjv(j/y) /Vij is the Rabi frequency 
for the A-photon resonance, T>ij = AEij/h — Aw is the 
detuning, fiij = j) is the dipole transition matrix 
element, A-Ejj = \Ej — E^] is the level spacing between 
the states i and j, Jjv is the Bessel function of integer 
order A, and y^j = dij-eo/fko with djj = fijj — flu- Note 
that Vij , , and T>ij depend on the number of photons 
A. In weak fields characterized by f2jj <C u> only tran- 
sitions with small A are important. Since the leakage 
is calculated from the maximum transition probabilities 
which are independent of pulse shape in weak fields jl^l 
the optimized WPs obtained using rectangular pulse are 
also valid for other pulse shapes. 

For the microwave-driven gate considered here, the mi- 
crowave is resonant with a pair of computational states. 
According to Eq. Q , the leakage is suppressed if Vij 3> 
flij holds for all unintended transitions. This condition 
may be satisfied by using sufficiently weak fields and/or 
qubits with proper spectroscopic properties. However, 
the use of exceedingly weak fields will make the gate very 
slow. Even so the unintended transitions will still cause 
large leakage at or near resonance where T>ij ~ and Vij 
~ I regardless of intensity of the fields. Thus the leak- 
age can be greatly reduced by setting the level spacings 
of the qubits substantially detuned from the microwave 
frequency and/or the dipole transition matrix elements 
sufficiently small for all unintended transitions. 

To investigate the leakage in a multi-level qubit, we 
study what happens when the microwave acts on each 
of the computational states, which we call a component 
of the gate. For an n-bit gate, there are M — 2" com- 



ponents. The leakage of the «th component is defined 
as Vi = S/c-^fc' where Pj& is the maximum probability 
of occupying an undesired state k through all possible 
multi-photon transitions and the sum is over all the un- 
desired states including all non-computational states as 
well as those computational states to which no transition 
is intended. The leakage of the n-bit gate, rj, is then de- 
fined as rj = max(?7i, rj 2 , • • • ?7m) which is a function of the 
v WPs q\, q 2 , ■ ■ ■, and q v . The set of WPs, g° pt , , • • •, 
and q° pt is obtained by minimizing the leakage of gate, 
namely, r? (g° pt , q^S • • • , <?° pt ) = min [77 (ft, q 2 , ■ ■ ■ , q v )\. 

As an example, we apply ITA to optimize the WPs of 
coupled rf SQUID flux qubits for CNOT gate. Each rf 
SQUID consists of a superconducting loop of inductance 
L interrupted by a Josephson tunnel junction character- 
ized by its critical current I c and shunt capacitance C 
[24l | . A flux-biased rf SQUID with total magnetic flux $ 
enclosed in the loop is analogous to a "flux" particle of 
mass m — C<I>q, where <J>o = h/2e is the flux quantum. 
The Hamiltonian is h (x) — p 2 /2m+mw 2 ^ c (x — x e ) 2 /2— 
Ej cos (27rx). Here, x — <&/<I>o is the canonical coordinate 
of the "flux" particle, p = —ihd/dx is the canonical mo- 
mentum conjugate to x, Ej — hl c /2e = mw| /C /3L/47r 2 
is the Josephson coupling energy, (3l — 2ttLI c /& is the 
potential shape parameter, lolc — l/^LC is the charac- 
teristic frequency of the SQUID, and x e = $ e /§o is the 
normalized external flux. 

The coupled rf SQUID qubits comprise two rf SQUID 
qubits: a control qubit and a t arg et qubit coupled via 
their mutual inductance M 0, Haj- For simplicity, we 
assume that the two SQUIDs are identical: d = C, 
Li = L, and I C i = I c for i = 1 and 2. The Hamil- 
tonian of the coupled SQUID qubits is H(xi,X2) — 
h (xi) + h (X2) + h\2 (xi, x%), where Xi, x e i, and h(xi) 
(i = 1 and 2) are the canonical coordinate, normalized 
external flux, and Hamiltonian of the ith single SQUID 
qubit and h\2 is the interaction between the qubits given 
by h 12 (xi,x 2 ) = mu} 2 LC n(xi - x el ) (x 2 - x e2 ). Here, 
k = M/L is the coupling constant. Note that each 
SQUID qubit is a multi-level system. The eigenstate 
|n) and eigenenergy E n of the coupled qubits are com- 
puted by numerically solving the eigenvalue equation of 
H(x\,x 2 ) using the two-dimensional Fourier- grid Hamil- 
tonian method [ifj. They are functions of the WPs x e \, 
Xf.,2 , and k for given device parameters L, C, and I c 
[27I Eif . For weak coupling k <^ 1 , the eigenstate of 
the coupled qubits, denoted by \n) — can be well 
approximated by the product of the control qubit's state 
\i) and the target qubit's state \j), \n) = \i) \ j). When 
biased at x e i, x e2 « 1/2 t he p otential of coupled SQUID 
qubits have four wells 0, l25| . The lowest eigenstates in 
each of the four wells, denoted as |1) = |00), |2) = |01), 
|3) = 1 10) , and |4) = |11), are used as the computational 
states of the coupled SQUID qubits. For SQUIDs with 
L = 100 pH, C = 40 fF, and L = 1.2, the energy lev- 
els and level spacings versus x e2 and k are plotted in 
Fig. 1(a) — (d), respectively. It is shown that both the 
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FIG. 1: Energy levels and level spacings of the coupled rf 
SQUID flux qubits. (a) the energy levels and (b) the level 
spacings vs. x e 2 for x e \ = 0.499 and k = 5 x 10 -4 , (c) the 
energy levels and (d) the level spacings vs. k for x e \ — 0.499 
and x e 2 = 0.49985. The lines 1 to 6 in (b) and (d) repre- 
sent the level spacings AE12, A.E13, AE14, AE23, AE24, and 
AE34, respectively. 



energy levels and level spacings are sensitive to the WPs 
x e2 and k. When x e2 < 0.49877 or x e2 > 0.50123 in 
Fig. 1(a) or k > 7.5 x 10~ 3 in Fig. 1(c) the energy level 
structures become quite complicated. Fig. 1(b) and Fig. 
1(d) also show that in certain regions of the WP space 
the level spacings become crowded, meaning that they 
are degenerate or nearly degenerate, which may result in 
significant intrinsic gate errors. 

Controlled two-bit gates can be realized by ap- 
plying a resonant microwave pulse to the target 
qubit. The interaction between the microwave and 
the coupled qubits can be written as V (xi, x 2 , t) — 
muj 2 LC {(x 2 - x e2 ) + k(xi- x e i) + x m /2] x m , where 
x m {t) — x m o cos (ujt) is the magnetic flux (normalized to 
$0) coupled to the target qubit from the microwave with 
amplitude x m0 and frequency to. For the CNOT gate, 
a 7r-pulse with w — AE^/h, where AE 34 — \E± — E 3 \, 
is used. Populations of the states 1 10) and |11) are 
exchanged after the 7r-pulse only if the initial state is 
1 10) or 1 11) or a linear combination of them. Note that 
in this case in addition to the non-computational states, 
the undesired states also include the computational 
states 1 00) and |01). Using Eq. Q and following the 
steps described above, the leakage of the CNOT gate 
i] = max(r?|oo>,J?|oi),??|io))^|ii)) is calculated. In the 
calculations, up to 3-photon processes are included. The 
TV-photon processes for N > 3 have negligible effects. 
For the SQUID qubits with the previously given DPs, 
the leakage is a function of WPs x e ±, x e2 , and k. In Fig. 
2(a) we plot the leakage of the CNOT gate versus x e2 
and k for x e \ — 0.499 and x m o = 2 x 10~ 4 . It is shown 
that the leakage of the CNOT gate is much smaller 
when the coupled SQUID qubits are operated around 
x e2 = 0.4997 and k = 7.5 x 10~ 4 for x el = 0.499. 

To evaluate the results obtained from ITA, we 
performed an accurate dynamic calculation by solv- 
ing the TDSE of the coupled qubits idc n (r) jdr = 



FIG. 2: Leakage of the CNOT gate for x e i = 0.499 and x m0 = 
2 x 10 -4 . (a) the leakage vs. x e 2 and k obtained from ITA, (b) 
the leakage vs. x e 2 and k obtained from DM, (c) comparison 
of the leakage vs. x e 2 at k = 5 x 10 -4 , and (d) comparison of 
the leakage vs. k at x e 2 = 0.4993. 

J2 n ' Hnn' ( T ) c «' ( T ) f° r the probability amplitudes c n 
of the first 20 eigenstates \n) using the split-operator 
method |2^|, where r = wlc* and H^ n , = 
[E n 6 nn > + (n \V\ n')] /hiu L c E^. The probability of be- 
ing in the state \n) is |c„| 2 , from which the maximum 
probabilities on the undesired states and the leakage 
77 (x e i, x e2 , k) of the CNOT gate are calculated. The re- 
sults are shown in Fig. 2(b) with the color scale identical 
to Fig. 2(a). To have a more quantitative comparison, 
we also plot, in Fig. 2(c) and (d), two line-cut figures at 
places with the richest structures in Fig. 2(a) and (b). 
The good agreement between the results of ITA and DM 
demonstrates the validity of using ITA to minimize the 
CNOT gate leakage. 

One of the advantages of ITA is that it provides clear 
physical intuition and insight into the origin of intrin- 
sic gate errors and thus how to reduce the problem by 
selecting appropriate WPs. Furthermore, ITA is orders 
of magnitude faster than DM. For instance, the result 
shown in Fig. 2(a) took less than 3 hours to compute on 
a dual-2.8 GHz processor Dell PRECISION 650 worksta- 
tion while that shown in Fig. 2(b) took more than 300 
hours on the same computer. Denoting t$ and tt re- 
spectively the time needed to calculate the spectroscopic 
properties of an n-bit gate and that used to compute the 
probability evolution in DM for each component. The 
total times needed to compute the leakage by ITA and 
DM are approximately 77 ~ Tg and rrj ~ Tg + 2 n r^, 
respectively, where the factor 2" is the number of com- 
ponents for the n-bit gate. Hence, Tg < tt and the 
ratio td /ti « 2™£ increases exponentially with n, where 
£ = rx/rg. For weak fields, £ ^> 1 and the ratio td/tj 
is very large. For example, for the CNOT gate above 
C ~ 25 and tjj/t] ~ 100. Thus as the number of qubits 
increases the optimization of WPs of multi-bit quantum 
gates could become an intractable problem using DM. 

To demonstrate quantitatively the effect of different 
WPs on intrinsic gate errors, we plot in Fig. 3(a) and 
(b) the population evolution of the computational states 
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FIG. 3: Population evolution of the computational states 1 10) 
and 1 11) for the CNOT gates with different WPs: (a) x e2 = 
0.49985 and k = 5x 10~ 4 corresponding to the point "A" in 
Fig. 2(a) [also (b)], and (b) x e2 = 0.49897 and k = 5 x 10~ 4 
corresponding to the point "B" in Fig. 2(a) [also (b)]. In 
both cases, x m o — 2 x 10~ 4 and x e i = 0.499. 

1 10) and 1 11) for the CNOT gate operated with two dif- 
ferent sets of WPs corresponding to points " A" and " B" 
marked in Fig. 2(a) [also (b)], respectively. It is clearly 
shown that inversion from the initial state 1 10) (|11)) to 
the final state |11) (|10)) after the 7r-pulse is almost com- 
plete when operated at the point "A" but significantly 
incomplete at the point "B". We also computed the pop- 
ulation evolution of the coupled qubits with the initial 
states 1 00) and |01) using the same microwave pulse. The 
system essentially stayed in the initial states for the gate 
operated at the point "A" but leaked significantly to the 
non-computational states for the gate at the point "B". 
The quality of a gate can be described by gate fidelity 
F = Trace [pppi], where pp and pi are the physical and 
ideal density matrices after gate operation and the over- 
line denotes averaging over all possible initial states [13 • 
For the CNOT gates operated at the point "A" we ob- 
tain F A = 0.9997 which is very close to the ideal CNOT 
gate. In contrast, the gate operated at the point "B" 
has Fg = 0.8061 which is too large to be tolerated. Fur- 
thermore, at the point "A" the gate is about a factor of 
three faster than that at the point " B" . 

In summary, a very efficient method is proposed to 
optimize the WPs of microwave-driven quantum logical 
gates in multi-level qubits based on ITA. In this method, 
the leakage is estimated accurately from the spectro- 
scopic properties of the qubits and minimized by choos- 
ing appropriate WPs. This method is exemplified by 
optimizing the WPs of coupled rf SQUID flux qubits for 
minimal leakage of the CNOT gate. The result is in good 
agreement with that obtained from dynamic calculations. 
Compared to the conventional dynamic method the ITA 
not only provides physical insight into the origin of gate 
leakage but also reduces the computational time by more 
than two orders of magnitude for the CNOT gate. Fur- 
thermore, since the ratio of the speed of ITA to that of 
DM scales exponentially as 2" the ITA is scalable as the 



number of qubit n increases. Our calculation also shows 
that high intrinsic fidelity CNOT gate can be achieved us- 
ing microwave-driven rf SQUID qubits with properly se- 
lected WPs. Although the proposed ITA is only applied, 
as an example, to the rf SQUID qubits, it is also valid 
for other microwave-driven multi- level qubits. Therefore, 
ITA provides a much needed solution for the optimization 
of WPs of multi-bit gates implemented with multi-level 
physical qubits, for which DM is extremely time consum- 
ing or could even become intractable. 
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